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Abstract. With combined hopping parameter and strong coupling expansions, we calculate 
a dimensionally reduced Polyakov-loop effective theory valid for heavy quarks at nonzero 
temperature and arbitrary chemical potential. We numerically compute the critical endpoint 
of the deconfinement transition as a function of quark masses and number of flavours. We also 
investigate the applicability of the model to the low-T and high density region, specifically in 
terms of baryon condensation phenomena. 



1. Introduction 

Despite its brilliant achievements in tackling perturbative problems, analytical investigations of 
QCD can not do much when it comes to determining the structure of its phase diagram in the 
(T, fj,) plane, with T temperature and fi quark chemical potential. In nonpertubative regimes, 
the answers numerical lattice QCD calculations can provide are mainly limited to ~ 0: the 
notoriuos sign problem renders the traditional Monte Carlo sampling approach meaningless 
and one has to rely on alternative methods, whose efficiency typically degrades about ~ T: 
examples are using an imaginary chemical potential, performing a Taylor expansion around 
fi = and reweighting of ensembles generated at zero chemical potential. 

Yet, the relevant physics for heavy- ion collisions, as well as for compact-star astrophysics, 
takes place at far-from-zero fi: it would then be desirable to devise strategies able somehow 
to circumvent the problems. In this category fall the many effective approaches developed 
throughout the years, such as sigma models, (p)NJL models and so on. 

In particular, in characterising the T — > region of the (T, phase diagram, the phenomenon 
of Silver Blaze is expected based on very general considerations: there, the partition function 
- hence all interesting observables - have to be independent of the chemical potential, up to 
a critical value associated to the lightest baryonic state in the theory. In the path-integral 
representation of QCD, the Dirac eigenvalues display an explicit dependence on fj., therefore 
the Silver Blaze phenomenon requires a highly non-trivial cancellation between them. Due to 
the shortcomings of the numerical methods mentioned above, the phenomenon has not been 
reproduced successfully in lattice QCD so far; also on the analytical side, an explicit proof exists 
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only for the choice of an isospin chemical potential = ^u = —fJ'd, with a critical onset located 
at /i/ = [1]. 

It is perhaps possible to avoid the sign problem altogether, thus exploring numerically the 
cold dense region /x ^ T of QCD, by employing the so-called Langevin dynamics [2j: recently, 
numerical results for a scalar complex field have been obtained [3] and later reproduced in the 
context of a flux representation for the same theory [Ij; both techniques, however, are at the 
moment hardly applicable to QCD, and it is not at all guaranteed that they will ever be. 

In this work we present an effective, dimensionally-reduced Polyakov-loop model obtained 
by applying systematically strong-coupling and hopping-parameter expansions to the original 
QCD lattice formulation, valid as long as quarks are not too light and the gauge coupling /3 
is not too large; after briefly describing the theory (Section [2]), we apply it to the cold dense 
region (Section [2D and show that the signature of the Silver Blaze phenomenon can be read off 
the numerical results. A more detailed presentation of the theory can be found in [5l^; for the 
Silver Blaze property in this context, see also [7]. 



2. The effective theory 

The model lives on a three-dimensional lattice, with complex scalars as per-site degrees of 
freedom, representing the traced Polyakov loops = Tr {Wx),Wx S SU{3). The associated 
effective action consists of a pure-gauge term and a fermionic contribution; the former comes from 
application of strong-coupling methods to the original four-dimensional Yang- Mills gauge action, 
and the latter is obtained with a hopping-parameter expansion, starting from the standard 
Wilson expression for the fermionic term in lattice QCD. In general, then, we have 

Zesi- ..) = j (Ild^.e^^^^)) , / dLe^(^) = f^O p\<t> e^^^"^^''^^) , (1) 

v ' 

[DL] 

where the parametrisation of the scalar L is made explicit (see [5]). The partition function 
Z depends on effective couplings which are in turn some functions of the original couplings 
appearing in the four-dimensional lattice QCD action. 

Let us now sketch how the pure-gauge part of the effective action is obtained [5] starting 
from the four-dimensional Euclidean Yang-Mills partition function at finite temperature (with 
time extent A',- lattice spacings) 

Zsu(,z),M = j m,] exp {A ^(Tr^7p + T:[/t)} , /3 = ^ = l8u + . . . . (2) 

By integrating out the spatial degrees of freedom and subsequently applying a strong-coupling 
expansion, one gets a tower of interactions in the effective theory, the leading one of which is 
simply a "spin-spin" term connecting nearest neighbours in the fundamental representation (the 
superscript "s" denotes the compliance of the terms to the centre symmetry requirement): 

- S',s = XiSf + \2SI • • • ; Ai^i^ = Ai(^, Nr)Y.{L,L* + L*Lj) ; Ai>i = o(Ai) for /3 ^ . (3) 

(ij) 

The precise form of the map Ai(/3,A^t-) comes from an order-by-order enumeration of strong- 
coupliiig graphs having the two Polyakov lines as boundary: we calculated the series up to order 
u^^ [5jo Moreover, one can resum higher powers of (certain classes of) graphs and improve the 
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effective theory to a "nonlinear" form 

Zeff = / [VL] + 2AiReL,L*) ; (4) 

the important point is that this model reproduces the deconfinement transition of the original 
theory, Eq. ([2]), and once the critical point (that still falls within the range of applicability of the 
strong-coupling methods) is known, it can be translated back to a table /3c(iV^) with a sufficient 
precision and in a wide enough range of to allow for a sensible continuum extrapolation of 
the physical deconfinement point Tc |5]. 

The next step is the inclusion of fermions in the theory: we perform an expansion in 
the hopping parameter k, therefore the masses (we consider A'^^ degenerate quarks of mass 
M = (1 — 8K)/{2aK)) must be large enough to guarantee convergence. The Wilson action for 
quarks is then written as 

oo i 

-S, = -NfJ2T^iH[UY) (5) 

with K is the hopping matrix. Its structure is such that the non-zero contributions after 
integration of the gauge fields are given by various kinds of closed loops of length l\ each of 
their links carries a factor k (and, if the chemical potential is turned on, an additional factor 
g±«M in the temporal direction): the expression above is then also an expansion in powers of k. 
The general form of the fermion contribution to the effective action is analogous to an external 
field in a spin system (here a denotes the explicit breaking of centre symmetry) H 

oo 

= -2A^/^(/ii(u,K,/i,iV05f +7Ii(^x,A^,/i,iV.)5,"'^) ; 7^,(-/i) = /ii(+//) , (6) 

i=l 

with the i-order refiecting increasingly subdominant terms. A closer inspection of the leading 
term SI reveals that, once again, a partial resummation of higher powers of the same graphs is 
possible: one finds that 

-SSj = log[ndet(l + /iiW^,)'^^det(l + 7^iT^j)'^^] +••• ; /ii = (2Ke'^'^)^- (1 + • • •) , (7) 

X 

The expression for h\ received also contributions from various types of higher-order graphs; 
including these, and the first subleading fermionic term S'f, one can write the model as [6l[7]: 



eff = liPL] n (l + 2AReL,L*)ndet[(l + /iiH^,)(l + /iiTyJ) 



IN, 



(we rename Ai — )• A and ignore higher-order corrections). Despite the appearance of Wx and 
in the above, can still be expressed, in practice, in terms only of the 0, (/> parameterising 
Lx in Eq. ([1]). The effective couplings are given by: 

/ii = cfl + eK^iV, "'^^ +■■■) ; (9) 
V \ — u ^ 

9 N 

= C^y(l + 2^--) ; (10) 
C = (2^6'^^)^- . (11) 



^ The next term seems an ordinary nearest-neighbour interaction, but is actually non-centre-symmetric, Eq. ((Sjl. 
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Figure 1. Left: Nf = 2 phase diagram in the space of temperature, chemical potential and 
meson mass (here called M^). Tq is the pure- gauge transition point; the blue line is the critical 
line. Right: T-^ phase diagram for the choice Nj = 2, M^d/T ~ 8.68; the vertical axis is 
normalised to the ^ = transition temperature for this mass choice. 



The meson and baryon masses are also evaluated: 

aruM = -21n(2K) - - 2Ak^— h . . . ; aniB = -31n(2K) - ISk^ — V ... ■ (12) 

1 — n 1 — n 

The model, as formulated in Eq. ([8]), is well suited to different simulational strategies: 
we could confirm that the standard Metropolis approach, a flux-representation-based worm 
algorithm, and a complex-Langevin implementation all agree with each other: however, one has 
to keep in mind that the first has to rely on a reweighting procedure for |U > (which does not 
hinder its scope significantly for not too large spatial volumes) and the second ceases to be of use 
if the /i2 term is included, while the third does not have a sign problem at all and is thus the best 
choice to deal with large-volume, finite-// ensembles. All of these algorithms, when compared to 
traditional lattice QCD approaches, are far less expensive in terms of required system resources 
and CPU time. 

An important application of this effective approach concerns the mapping of the critical 
surface in the (M^^, Mg, fi/T) space. Due to the heavy-quark validity of the model, the second- 
order surface that could be located is the one associated to the upper-right corner of the Columbia 
plot: there, a nice agreement between the findings at zero, real positive and imaginary fj, was 
verified and the correct universality class of the transition was confirmed. Moreover, the result 
at zero chemical potential for A'^^ = 1,2,3 closely reproduces those found in ordinary QCD 
simulations e.g. in [21 E]. The shape of the critical surface at all chemical potentials can be 
summarised in the parametrisation [B]: 

(A^//.ie-/^).. = '-^^ ; - e-^^V2 , M, ^ oo ; (13) 

the phase diagram for heavy-quark QCD can then be determined in (y, ^f-, ^), and, by taking 
slices of it, a "heavy fermion" version of the familiar /i-T phase diagram can be drawn (Fig. [1]). 

3. Low-temperature regime 

We now turn to investigate the cold and dense regime of QCD with this effective theory. In 
the (analytically solvable) static strong-coupling limit, one finds that the Silver Blaze property 
holds as T ^ 0, with the quark number density approaching a step function centred at 
3/i = niB = — 31og(2K) (the last equality is exact in this limit): 

a\'^2N,&{fiB-mB) ; (14) 
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Figure 2. Continuum-extrapolated baryon density, in units of function of the baryon 

chemical potential (also measured in units of ms) for various temperatures. 



one sees then also that a saturation phenomenon is expected (this is connected to the 
resummations performed in deriving the form of Zcg). The next step is to perform numerical 
computations with the full expression Eq. ([8]). Taking advantage of the fact that Nr is now 
merely a parameter in the effective couplings' maps, we can tune it to very large values in the 
hundreds: this effectively realises A — > thus simplifying the model; by keeping k 10^'^, 
moreover, we can still trust the pure- Yang-Mills scale-setting prescription [1^. In this way, we 
can measure numerically the baryon density as a function of the baryon chemical potential with 
m-B = 30 GeV and four different temperatures from 20 MeV to 2.5 MeV [7]. 

The baryon density ub can be expressed as a meaningful physical quantity after a careful 
continuum limit is taken: we could generate data at nine different values of the lattice spacings 
a and perform a fit to the continuum limit; Uf, is measured as 



The results are shown in Fig. [2l the step-like shape of the density is less and less smoothed- 
out as r ^ (and anyway the //-scale is extremely tiny); besides, as the baryonic chemical 
potential hits the baryon mass the expected saturation to a common, roughly T-independent 
density is seen, and this value, when expressed in appropriate units of m^, is less than a factor 
two off the real-world nuclear matter, despite this model being limited to extremely massive 
quarks. This seems to suggest that the core features of nuclear condensation might be largely 
independent of whether the fundamental constituents of the baryon are heavy or light. 
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